Relaxation, dephasing, and quantum control of electron spins in double quantum dots 
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Recent experiments have demonstrated quantum manipulation of two-electron spin states in dou- 
ble quantum dots using electrically controlled exchange interactions. Here, we present a detailed 
theory for electron spin dynamics in two-electron double dot systems that was used to guide those 
experiments and analyze the results. Specifically, we analyze both spin and charge relaxation and 
dephasing mechanisms that are relevant to experiments, and discuss practical approaches for quan- 
tum control of two-electron systems. We show that both charge and spin dephasing play important 
roles in the dynamics of the two-spin system, but neither represents a fundamental limit for electrical 
control of spin degrees of freedom in semiconductor quantum bits. 
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Electron spins in quantum dots represent a promis- 
ing system for studying mesoscopic physics, developing 
elements for spintronicsi^, and creating building blocks 
for quantum information processing 3 .! 4 -^. In the field 
of quantum information, confined electron spins have 
been suggested as a potential realization of a quantum 
bit, due to their potential for long coherence times 7 -^ 9 -. 
However, the deleterious effects of hyperfine coupling 
to lattice nuclear spin a 10 : 11 ! 12 : 13 ' 14 : 15 ! 16 ! 17 ! 18 , as found in 
experiment a 19 ! 20 ' 21 ! 22 ' 23 , can severely limit the phase co- 
herence of electron spins. Thus, it is important to un- 
derstand dynamics of electron spin coupled to nuclei and 
to develop corresponding quantum control techniques to 
mitigate this coupling. 

Recent experiments by our group explored coherent 
spin manipulation of electron spins to observe and sup- 
press the hyperfine interactio n 20 ! 21 : 24 . In this paper, we 
present a detailed theory describing coherent properties 
of coupled electrons in double quantum dots that was 
used to guide those experiments and analyze the results. 
The theory includes hyperfine interactions, external mag- 
netic field, exchange terms, and charge interactions. 

Our approach relies upon an approximation based on 
the separation of time scales between electron spin dy- 
namics and nuclear spin dynamics. In particular, the 
time scales governing nuclear spin evolution are slower 
than most relevant electron spin processes. This al- 
lows us to treat the nuclear environment using a type 
of adiabatic approximation, the quasi-static approxima- 
tion (QSA) 11 : 16 . In this model, the nuclear configuration 
is fixed over electron spin precession times, but changes 
randomly on the time scale over which data points in an 
experiment might be averaged (current experiments ac- 
quire a single data point on ~ 100 ms timescales). We 
also consider the first corrections to this approximation, 
where experimentally relevant. 

In what follows, we start by reviewing the theory of hy- 
perfine interactions in single and double quantum dots, 
focusing on electrostatic control of electron spin-electron 



spin interactions. We then consider the role of charge de- 
phasing and charge-based decay in experiments involving 
so-called spin blockade, in which a simultaneous spin flip 
and charge transition is required for electrons to tunnel 
from one dot to another—. Consistent with the experi- 
ments, we find that blockade is reduced near zero mag- 
netic field over a range set by the average magnitude of 
the random Overhauser (nuclear) field. We then consider 
the effect of fast control of the local electrostatic poten- 
tials of double quantum dots, and show how this may be 
used to perform exchange gate a 3 ' 10 : 25 , and to prepare and 
measure two-spin entangled state a 21 ' 26 . Various limita- 
tions to the preparation, manipulation, and measurement 
techniques, due to nuclear spins, phonons, and classical 
noise sources, are considered. 

Theories that explicitly include quantum mechanical 
state and evolution effects of the nuclear spins both 
within and beyond the QSA have been considered by sev- 
eral authors (Refs.0l3lllIIillllllli). Dephasing, 
decoherence, and gating error in double quantum dots 
have also been investigated previousl y 27 ' 30 : 31 ; the present 
work develops the theory behind quantum control tech- 
niques used in experiments, connecting the previous gen- 
eral theoretical treatments to specific experimental obser- 
vations. The paper is organized as follows. Interactions 
of a single electron in a single quantum dot, including 
hyperfine terms, are reviewed in section [T] The quasi- 
static regime is defined and investigated, and dephasing 
of electron spins by hyperfine interactions in the quasi- 
static regime is detailed. This provides a basis for ex- 
tending the results to double quantum dot systems. We 
then develop a theory describing the two-electron spin 
states of a double quantum dot including the response of 
the system to changes in external gate voltages, and the 
role of inelastic charge transition a 32 ' 33 ' 34 . This is com- 
bined with the theory of spin interactions in a single dot 
to produce a theory describing the dynamics of the low 
energy states, including spin terms, of the double dot sys- 
tem in two experimental regimes. One is near the charge 
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transition between the two dots, where the charge state 
of two electrons in one dot is nearly degenerate with the 
state with one electron in each dot. The other is the far- 
detuned regime, where the two dots are balanced such 
that the states with two electrons in either dot are much 
higher in energy. 

In the remaining sections, we investigate situations re- 
lated to the experiments. First we consider spin block- 
ade near the charge transition, as investigated in Ref . [2(| 
considering effects due to difference in dot sizes and ex- 
panding upon several earlier, informal ideas. We then 
analyze approaches to probing dephasing and exchange 
interactions, showing how errors effect fast gate con- 
trol approaches for preparation and measurement of two- 
electron spin states, as well as controlled exchange inter- 
actions and probing of nuclear-spin-related dephasing, as 
investigated in Ref. l2l|. Finally, we consider limitations 
to exchange gates and quantum memory of logical qubits 
encoded in double dot systems^S^. 

I. HYPERFINE INTERACTIONS IN A DOUBLE 
QUANTUM DOT: A REVIEW 

We begin by reviewing the basic physics of hypcrfinc 
interactions for electron spins in single GaAs quantum 
dots. This section reviews the established theory for sin- 
gle quantum dot o 11 ' 12 ' 16 ' 35 and considers dynamical cor- 
rections to the model of Refs. [I0GJ. This model will 
be used in subsequent sections for the double-dot case. 
Additional terms, such as spin-orbit coupling, are ne- 
glected. Theory 7 -^ and experimen t 20 ' 36 ' 37 have demon- 
strated that spin-orbit related terms lead to dephasing 
and relaxation on time scales of milliseconds, whereas we 
will focus on interaction times on the order of nanosec- 
onds to microseconds. 

A. Electron spin Hamiltonian for a single quantum 

dot 

The Hamiltonian for the Kramer's doublet of the 
ground orbital state of the quantum dot (denoted by the 

spin-1/2 vector S) including hyperfine contact interac- 
tions with lattice nuclei (spins 7^ ,J ) i a 11 ' 38 

h = h le B cxt .s + h le Y, h<*j/s ■ (!) 

0,3 

where 7 e = g* fis/fr is the gyromagnetic ratio for electron 

spin S; sums are over nuclear species (/3) and unit cells 
(J). Correspondingly, bp is the effective hypcrfinc field 
due to species (5 within a unit cell, with 675 As = —1.84 T, 
&69 Ga = -1-52 T, and &7i Ga = -1.95 T for GaAs 3 ^, The 
coefficient ctj_p — wo|V'(' ? j,/3)| 2 is the probability of the 
electron being at unit cell j (with nuclear spin species 
vq is the volume of the unit cell (2 nuclei), and tp(r) 
is the envelope wavefunction of the localized electron. 




FIG. I: (a) A schematic potential and energy level diagram 
for a single quantum dot in which one electron is confined 
to the low energy spectrum of a three dimensional potential. 
Only the ground and first-excited states, each a Kramer's 
doublet, are shown, (b) The lowest orbital state has a spin- 
1/2 electron interacting with the lattice nuclear spins, (c) 
Effective magnetic field due to both external field and the 
nuclear field. When the external field is large, the transverse 
components of the nuclear field are neglected in a rotating 
wave approximation. 



It is convenient to rewrite the Hamiltonian using 
a collective operator for the nuclear spins, i? n uc = 

^pTlij a j,pl • This operator allows us to write the 
Hamiltonian as an electron spin interacting with an ex- 
ternal magnetic field, -B ox t, and an intrinsic field, B nuc : 



H eS = h le (B cxt + B nuc ) ■ S . (2) 



Several characteristic value o 11 ' 38 for this interaction 
are noted in Table |TJ The maximum nuclear field 
value (all spins fully polarized with value I = 3/2) is 
ho = Y^g bp{xplP) ^2 k where we have separated out 
the relative population of nuclear species, £75 As = f, 
^BfGa = 0.6, and X7i Ga = 0.4 for GaAs, removing the 
f3 dependence from the a>k,/3- This gives bo — 5.3 T. 
Second, when the nuclear spins may be described by a 
density matrix p — 1/ (21 + l) N (infinite temperature ap- 
proximation for N nuclei), the root-mean-square (rms) 
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strength of the field 1 is 



decay as: 



'<|5 nuc H/3= E&K,^l f/3 ' fc l 2 )/ 3 ( 3 ) 
V 0,k 

£>& 2 )I(I + 1)V3 J 'd*r\i>{r)\± (4) 

(5) 





= hi/VN 

where we have replaced J2j v o with J d 3 r. 
The characteristic strength parameter is hi = 

^21(1 +l)/3J2p x P b } = 4 -° T for GaAs ' and N 
is defined as the number of nuclei with which the elec- 
tron has significant overlap, i.e., N — 2/[J d 3 r \i/j(r)\ 4 vo]. 
These numerical values are specific for GaAs quantum 
dots. Dots in other materials with non-zero nuclear spin 
may be described by similar parameters: a maximum 
field strength parameter ho and a rms field strength 
parameter, B nuc = h±/yN. 



B. The quasi-static approximation for nuclear 
spins 

By writing the Hamiltonian (Eq. \2§ with nuclei as an 
effective magnetic field, we have implicitly indicated that 
the field may be considered on a similar footing to the 
external magnetic field. In other words, the operator 

Bnuc may be replaced by a random, classical vector B nuc , 
and observables may be calculated by averaging over the 
distribution of classical values. The distribution in the 
large N limit is 



P{B) 



1 



(2nBl u 



)3/2 



cxp(-(B ■ B)/2B 



2 ) 

nuc/ 



(6) 



This is the quasi-static approximation (QSA) used in 
Refs.mil 2 : we assume that over time scales correspond- 
ing to electron spin evolution, the nuclear terms do not 
vary. 

In terms of dephasing, we cite the results of Ref. [Til . 
At zero external magnetic field in the Heisenberg picture, 

the electron spin S evolves to: 



(S(t)) t 



5(0) 



(l + 2(l-( 7e B nuc t) 2 )e-^ B - t ) 2 ). (7) 



On the other hand, at large external magnetic fields, S z 
is conserved, but transverse spin components (e.g., S x ) 



1 A different convention than that of Ref. [Tl] is used: gaussians 
are always described by their rms, rather than vl times their 
rms. Thus our value of B nuc = A.g/y/2, and similarly for other 
values to follow. 

2 The QSA is designated the quasi-stationary fluctuating fields 
approximation in Ref. Illl . 
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(8) 



A time-ensemble-averaged dephasing time due to nuclei 
in a single dot at large external magnetic field (e.g., dot 

i) is 



rp* |_ 



(9) 



This definition is appropriate when considering the decay 
of coherence of a single electron in a single quantum dot. 
Generalizing to all field values, for times longer than 



T. 



2 ) 



S(t > T 2 *) = ((S(t = 0) • n)n), 



(10) 



is the average electron spin value, averaged over a time 
r = 2n/u!. 

At low magnetic fields, the QSA is valid up to the single 
electron spin- nuclear spin interaction time (0(hN/A)), 
which is of order microseconds^ 1 . In contrast, at large 
external fields, the regime of validity for the QSA is ex- 
tended. Terms non-commuting with the Zeeman interac- 
tion may be eliminated (secular approximation or rotat- 
ing wave approximation), yielding an effective Hamilto- 
nian 



HcS — h r y e (B cxt + B^ UC )S Z 



(11) 



The z axis is set to be parallel to the external mag- 
netic field. Corrections to the QSA have a simple in- 
terpretation in the large field limit. As the Zeeman en- 
ergy suppresses spin-flip processes, we can create an ef- 
fective Hamiltonian expanded in powers of l/B cxt us- 
ing a Schrieffer-Wolf transformation. In the interac- 
tion picture, we write the corrections to H Q g by setting 
B* uc = B z + SB z (t), where B z is the QSA term, and 
8B z (t) are fluctuations beyond the QSA. When the num- 
ber of nuclei N is large and fluctuations small, we approx- 
imate 5B Z by its fourier transformed correlation function: 



(6B z (t + T)5B z (t)) = / duj S(ui)e 



(12) 



where S(lo) has a high frequency cutoff 7 <C 7e-B nuc . 
The form of S{uj) depends on the detailed parameters 
of the nuclear spin Hamiltonian and the nuclear spin- 
nuclear spin interactions, and in general requires a many- 
body treatment. A variety of approaches have been used 
to successfully estimate these correction o 17 ' 18 i 27 ' 28 i 29 ' 39 . 
Any approach with an expansion in inverse powers of 
the external field is compatible with our assumption of 
S(uj), provided that the number of nuclear spins is suf- 
ficiently large that Gaussian statistics may emerge. In 
contrast, the validity of the QSA in the low-field regime 
remains unproven, though recent simulations 4 ^ suggest it 
may break-down before h\^N j A timescales. 
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Type 


Time 


Energy 


Magnetic field 


Typical value 


Charge 










Charging energy 




E c 




5 meV 


Orbital level spacing 




hu>o 




1 meV 


Single dot two-electron exchange near _B ext = 




J 




300 fteV 


Double-dot tunnel coupling 




T c 




10 A*eV 


Double-dot inelastic tunneling 




r(e) 




0.01-100 neV 


Electron spin 










Larmor precession 


t L 




£>cxt 


0-200 MeV 


Fully polarized overhauser shift 




AI 


h 


130 fieV 


IK atifiAvn 1 (morninicov cViitt 
1 XXailUUIll 1 U V til llclLLOcI OllllL 


1 2 


h'v R 


f> — Jn_ 


n i-i /<pv 


Nuclear spin species f3 










Larmor precession 


tnL,f3 




Bcxt 


0-100 neV 


Knight shift 


tic,p 


fi-y e B nuc \0j sa 7 ^ ut 


7eS„uc 

1/3 -/N 


0.1-10 neV 


Dipole-dipole interaction (nearest neighbor) 


tdd 


vn 


v n 


0.01 neV 



TABLE I: Time, energy, and magnetic field scales for electron and nuclear spins in a single and double quantum dots, from 
fast to slow. 



C. Hyperfine interactions in a double quantum dot |(1, 1)5) can be written as 



We consider standard extensions to the single electron 
theoretical model to describe the case of two electrons 
in adjacent, coupled quantum dots, by considering only 
charge-related couplings, then including spin couplings. 
The relevant states are separated electron states, in which 
one electron is in each quantum dot, and doubly-occupied 
states, with two electrons in one of the two dots. 

The doubly-occupied states are assumed to be singlets 
(appropriate for small perpendicular magnetic field)^. 
The higher excited states that are doubly-occupied are 
triplets with a large energy gap J . This singlet-triplet 
energy gap for doubly-occupied states facilitates elimi- 
nation of the spin interactions and the doubly-occupied 
triplet states. Furthermore, by controlling the relative 
potential e of the two quantum dots using electrostatic 
potentials applied by external gates, the ground state can 
be changed from one of the doubly-occupied states to one 
of the separated electron states (far-detuned regime on 
the other side of the charge transition)^. Electrostatic 
control of the double dot Hamiltonian will be analyzed 
in more detail in sections ILT1 IIV1 



Formally, we eliminate all but one of the doubly- 
occupied states following the prescription of Ref. H3- We 
include the doubly-occupied state (0,2)S, where (n;,n r ) 
denotes number of electrons in left, right dots respec- 
tively, and S denotes a singlet of electron spin, in addition 
to the singlet and triplet manifolds of the (1,1) subspace. 
For notational convenience, we set e = to occur at the 
avoided crossing between (1,1) and (0,2) in Fig. [2] There 
is an avoided crossing at e = for the spin singlet man- 
ifold due to quantum mechanical tunneling T c between 
the two quantum dots, while the spin triplet manifold 
is unaffected. The Hamiltonian for the states | (0,2)5), 




(13) 



As the tunneling coefficient and external magnetic field 
are assumed constant, we will choose T c to be real by an 
appropriate choice of gauge 

For a slowly varying or time- independent Hamiltonian, 
the eigenstates of Eq. [13] are given by 

5^ = cos6»|5) + sin 6»| (0,2)5) (14) 
G^ = - sin6» |5) + cos 9 |(0, 2)5) . (15) 

We introduce the tilde states as the adiabatic states, with 
the higher energy state. The adiabatic angle is 9 — 
arctan(- 
are 



2T C 



), and the energies of the two states 



E § = 



Eg 



-ytan(0) 
T 

— tan(7r/2 — 



(16) 
(17) 



When 



5 



ITJ, 9 



0, the eigenstates become 



|5), G) -> |(0,2)5). For e » |T C |, 9 -> tt/2, 



5 -> |(0,2)5) 



and the eigenstates are switched, with 



and 



G 



1 5). As will be discussed later, controllably 

changing e allows for adiabatic passage between the near 
degenerate spin states |5) , \T m ) (far detuned regime) to 
past the charge transition, with | (0,2)5) as the ground 
state (e 3> \T C \). This adiabatic passage can be used for 
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FIG. 2: (a) Charge stability diagram for a double dot system. 
Double dot occupation is denoted by (ni,n r ). The detuning 
is parameterized by e, and the far-detuned regime (light blue) 
and charge transition (yellow) are shown, (b) Schematic of the 
double well potentials along one axis (a;) with tight confine- 
ment in the other two axes (i.e., y and z). In the far-detuned 
regime, the (1,1) charge states are the ground state, while in 
the charge transition regime, (0,2) can be the ground state. 
Triplet states are indicated in red, while electron charges are 
indicated in orange, (c) Energy level structure of the double 
dot system as a function of detuning. From left to right the 
lowest energy charge state as a function of e is (2,0), (1,1), and 
(0,2). The detuning at the middle of the graph corresponds 
to e = —E c /2, where E c is the charging energy of a single 
dot. The three (1,1) triplet states (shown in red) are split by 
Zeeman energy. 



singlet generation, singlet detection, and implementation 
of exchange gatesSi. 

We now add spin couplings to the double dot sys- 
tem, including both Zeeman interactions and hyperfme 
contact coupling. Two effective Hamiltonians, one for 
e <C — \T C \ and one for e ~ are developed. Our ap- 
proach is similar to that of Ref. |27|, and we include it 
here for completeness. The spin interactions in a double 



quantum dot for the states \T„ 
e< -\T e \ as 



S ) may be written for 



hf.tot 



l 

hf,eff 



H^ eS -J(e)\S)(S\ (18) 



where / and r refer to left and right dot, respectively, the 
nuclear fields are determined by the ground orbital state 
envelope wavefunctions of the single dot Hamiltonians 
(see Fig.©, and J(e) = -E§(e). 

Reordering terms simplifies the expression: 



Hi 



hf.tot 



hj e [B ■ (S l + S r ) + dB ■ (S 1 - S r )} 



J(e)\S)(S\ 
(19) 



(b) 



s> 



t- 



FIG. 3: A double quantum dot in the (1,1) configuration, (a) 
Schematic of the two-electron wavefunction in the far-detuned 
regime interacting with lattice nuclear spins, (b) Electron 
spins in the left and right dots interacting with their respective 
effective nuclear fields in the quasi-static approximation. 



(a) 
|T >_ 




FIG. 4: (a) Levels in the far-detuned regime including all 
couplings of Eqn. 1201 (b) Levels near the charge transition; 

the \T+) *-* is near resonance, with the coupling between 
|T+) and \S) indicated, as per Eqn. 1221 



i+S„ 



and differ- 



with an average field B — B cxt H j~ 

cncc field dB = (B nuc . } - B nuc , r )/2. The form of Eqn. EH 
indicates that terms with B and J(e) are diagonal in to- 
tal spin and spin projection along B, creating a natural 
set of singlet and triplet states. However, the term with 
dB breaks total spin symmetry, and couples the singlet 
to the triplet states. 

We can now write Eqn. [115] in matrix form in the basis 
{\T+),\T ),\T-),\S)}, 



( 



H = h 7e 



B z 









-dB, 






-B z 

-dB x -\-idB v 



idB u 



V2 

-dB z 



\ 



V2 

- a ** "vsT" 

The corresponding level structure is given in Fig.^i. We 
have implicitly assumed the QSA in writing this Hamil- 
tonian by defining the axis of spin up and down as B, 
which is a sum of the external field and the average nu- 
clear field. If the nuclear field fluctuates, those terms can 
contribute by coupling different triplet states together. 

With no external magnetic field, all states couple to the 
singlet, and solving the dynamics requires diagonalizing 
the 4-by-4 matrix of Eqn.[20l However, at finite magnetic 
field, a large Zeeman splitting (which sets B z » B nuc ) 
allows us to separate the system. The far-detuned regime 
only has transitions between the m s = states; in this 
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basis ({|Tq) , |5)}), the matrix becomes 




(21) 



This two-level system has appropriately straightforward 
dynamics, and we investigate it in some detail below. 

Near the charge transition and at finite magnetic field, 
another coupling can occur, this time between \T + ) and 

S ) . This resonance corresponds to the adiabatic singlet 



5y having an exchange energy J(e) close to the Zeeman- 

split triplet's Zeeman energy, E z = hrf e B z . We note that 
external magnetic field for GaAs will be negative in this 
context (due to the negative electron g- factor = -0.44). 

Written in the basis {|T+) , S)}, the Hamiltonian is 



ip— flop 



B z 

idB-u 



dB x —idB v 



V2 



■ COS ( 



\/2 

-J(e)/le 



We have subscripted Eqn. with "flip-flop'' 



(22) 
indicating 



5) and |T + ) result in the flipping of 



that flips between 
a nuclear spin, which can be seen by identifying dB + = 
dB x + idBy — (i?„uc,i,+ — Sjiuc,r,+)/2. 

Because the S-T + resonance leads to spin flips and 
eventual polarization of the nuclear field, the QSA will 
not be valid if appreciable change of field occurs, and 
the overall dynamics may go beyond the approximation. 
This has been examined experimentall y 22 ! 42 and theoret- 
icall y 43 ' 44 ! 46 for some specific cases. While the discussion 
to follow mentions this resonance, it will focus on the zero 
field mixing of Eqn. ll9l and the far-detuned regime's finite 
field mixing of Eqn.JHU We remark that Eqns.l2"Tlandl2"2l 
have been previously derived outside of the QSA 27 . 

We have now established that in the far-detuned 
regime, the relevant spin interactions are limited to dy- 
namics within the singlet-triplet subspace and deter- 
mined by the Hamiltonian in Eqn. 1211 Similarly, near 
the charge transition a resonance between |5) and |T+) 
may be observed; as this resonance allows for nuclear 
spin polarization, it may only be partially described by 
the QSA, and we do not consider its dynamics in de- 
tail. However, we note that in the absence of nuclear 
spin polarization the resonance occurs when the Zeeman 
splitting of the external field equals the exchange energy, 
J(e). Thus, if the Zeeman energy is known, measuring 
the position of the splitting gives a map between external 
parameters and the actual exchange energy. 



II. NUCLEAR-SPIN-MEDIATED RELAXATION 
IN DOUBLE DOTS 

In this section we consider the case in which the ground 
state of the system is |(0, 2)5) and the low-lying excited 



states are the (1,1) states, {|5) , \T m )}. This situation oc- 
curs in dc transport when the system is in the spin block- 
ade regime, where transitions from \T m ) to 1(0,2)5) are 
suppressed because they require both a spin and charge 
transition. Previous theoretical work for two electron 
systems has focused on triplet and singlet decay of two- 
electron states in a single quantum dot^£; in contrast, 
the present analysis deals with a double quantum dot 
system where the electrons can be well separated. Con- 
trary to more general spin blockade calculations 4 ^ and 
experiments, the present work is focused entirely on the 
rate limiting step of blockade: the spin flip followed by 
charge transition within the double quantum dot. Sev- 
eral group o 36 i 37 ' 49 have studied spin relaxation between 
Zeeman split spin states at high magnetic field (B >4 
T). The measured relaxation rates were found to scale as 
B , consistent with a spin-orbit mediated spin relaxation 
process^. Similarly, single dot measurements of triplet- 
singlet relaxation when J ^> 7 e -B n uc (i-e., when the effect 
of nuclei is small) indicate long lifetimes, likely limited by 
similar spin-orbit mediated mechanisms or cotunneling to 
the lead a 33 ' 50 . On the other hand, at low field and small 
exchange, when the splitting between spin states becomes 
comparable to B nuci the hyperfine interaction dramati- 
cally increases the spin relaxation rate. Recent experi- 
ments have measured spin relaxation between nearly de- 
generate singlet and triplet spin states in this regim o 20 ' 51 . 
Experimental techniques are discussed in Refi^ 2 -, and a 
full analysis of the field and energy dependence of the 
relaxation rate is discussed in Refi2£. We only briefly 
outline the salient features of experiment, focusing in- 
stead on developing a more rigorous basis for the theory 
of previous, published work. 

Experiments are performed near the two-electron 
regime with very weak tunnel coupling so that T c is 
slower than the pulse rise times (T c -C 1/xeV). Pulsed- 
gate techniques are used to change the charge occupancy 
from (0,1) to (1,1) to (0,2) and back to (0,1). In the (1,1) 
charge configuration with weak interdot tunnel coupling 
the |5) and \T m ) states are nearly degenerate. Shift- 
ing the gates from (0,1) to (1,1) creates a mixture of all 
four states, |5) , |I^i=-i,o,i) by loading an electron from 
a nearby Fermi sea. Then, the system is rapidly (non- 
adiabatic with respect to tunnel coupling T c ) shifted to 
the (0,2) regime, with 1(0,2)5") as the ground state. In 
this rapid shift procedure, the singlet |5) does not adia- 
batically follow to the doubly-occupied singlet | (0,2)5), 
but instead follows the Zener branch of the avoided cross- 
ing and stays in |5), as is illustrated in Fig. [5] 

Past the charge transition, when the adiabatic basis 

G ) is an appropriate representation of the system, 



S 

it is possible for the system to experience inelastic de- 
cay from the excited state G) to the ground state 5 
via charge coupling, e.g 

> 



for 



T, 



is E d 



Gj to the ground state 

to p honons. Th e energy gap 
= y/e 2 + 4|T C | 2 w e. Inelas- 



tic decay near a charge transition in a double quantum 
dot has been investigated in great detai l 32 ' 34 ' 53 , and we 
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FIG. 5: (a) Energy level structure as a function of detuning. 
Coupling to a Fermi sea with kbT 3> 7 e |-B C xt| leads to equal 
filling of all four low energy states in the far-detuned regime 
(labelled start). Then e is changed rapidly with respect to 
the tunnel coupling, leading to all four spin states still in the 
(1,1) charge configuration. The time spent waiting in this 
configuration results in slow decay of the metastable (1,1) 
states to the |(0, 2)S) state, (b) Nuclear spins couple between 
the eigenstates of exchange, and slow inelastic decay at a rate 
F proceeds from \S) to |(0, 2)S). (c) The same process, but for 
the eigenstates of the double dot Hamiltonian. The m a = 
states are equal superpositions of \S) and decay rapidly to 
|(0, 2)S), while the \m 3 \ — 1 states are only weakly mixed at 
large magnetic field with the \S), resulting in slow decay to 
1(0,2)5). 



do not seek to reproduce those results. Instead, we note 



that the decay from the excited state 
state 



Gj to the ground 

Sj is well described by a smoothly varying, energy- 
dependent decay rate T p h(e). Incoherent population of 
G^ by absorption of a thermal phonon is suppressed as 

long as 2T C > k^T, which is satisfied for T c = 0.01 meV 
and T = 100 mK. 

Finally, we can combine the coherent spin precession 
due to interaction with nuclear spins with the charge- 
based decay and dephasing mechanisms to investigate 



relaxation of \T m ) states to the state 



Of partic- 



ular interest is the regime past the charge transition, 

e ^> T c , where « |5) becomes nearly degenerate 

with \Tq), as was studied in the experiment of Ref. I20I . 
An effective five-level system is formed with the levels 
{|T+) , \T ) , |T_) , IS")} described by the spin Hamiltonian 
of Eqn.[2Dl while inelastic decay from |5) to |(0, 2)5) (the 
fifth level) is possible at a rate T(e), as shown in Fig. [SJd. 

To analyze this process, we start with the Louivillian 



superoperator that describes inelastic tunneling: 

p = i[p, H/h] +r(e)/2[|5X5|p + p|5X5| 
-2|(0,2)5)(5|p|5)((0,2)5| 

where 



(23) 



H = h le [B t -S l +B r -S r ] (hl) ~he\(0,2)S)((0 1 2)S\ . (24) 

I and r indicate left and right spins for the (1,1) charge 
space. Assuming the nuclear field is quasi-static (QSA), 
we can diagonalize H. The eigenstates are the ground 
state, |(0,2)5), and (1,1) states with spin aligned and 
anti-aligned with the local magnetic fields, B l ' r = B oyLt + 
£?nuc- We write these eigenstates as \s,s') = |s) ; <§5 \s) r , 
where s,s' — ±1/2 are the eigenvalues of the spin pro- 
jection on the fields of the I and r dots, respectively. The 
eigenvalue for | (0,2)5) is E G = —e and the other four 
eigenstates \s,s') have energy 



E s , s , =sj e \B l \ + s'^ e \B r \ 



(25) 



In considering decay from the energy eigenstates of 
the nuclear field, \s,s'} to |(0,2)5), we eliminate rapidly 
varying phase terms, e.g., 1 1/2, -1/2) (-1/2, 1/2|. This 
is appropriate provided that the inelastic decay mecha- 
nism, r(e), is slow in comparison to the electrons' Lar- 
mor precession in the nuclear field -B n uc- In this limit, 
each state |s, s') decays to |(0, 2)5) with a rate given by 
r(e)| (s, s' |(0, 2)5) | 2 , as indicated in Fig. [5b. A detailed 
analysis is given in appendix [A] For convenience, we 
write c s , s / = (s,s' \ (0,2)5). 

Starting with a mixed state of the (1,1) subspace (as in 
Ref. l2(il ). we can find analytical expressions for the time 
evolution of the density matrix of an initial form p(t = 
0) = J2 S s' I s ' S 'X S > s 'l/4- This initial state corresponds 
to a mixture of the four (1,1) spin states. The charge 
measurement distinguishes only between (1,1) states and 
|(0,2)5); accordingly, we evaluate the evolution of the 
projector for the (1,1) subspace Pn = J2 S s , \s, s')(s, s'\. 
In particular, 



-r(e)|c 4 



-l 2 */2- 



-r(e)|c-t 



-i 2 72 



(26) 



Finally, we must average over possible initial nuclear 
spin configurations to find the measured signal. This 
means evaluating (Pn(t)) n uc, a difficult task in general. 
However, 

(PiiWW w e -r( £ )<|c ++ | 2 )» uct/2 + e -r( £ )(| c+ _|^) nuct/2 _ 

In this approximation, we replace the average of the 
exponents with the average values for the coefficients 

|c ++ | 2 ,|c_| | 2 . The validity of this approximation can 

be checked with numerical integration, and for the range 
of parameters presented here the approximation holds to 
better than 1%. 
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FIG. 6: The product Iil r as a function of external mag- 
netic field B in units of average nuclear field, B nuc = 
(Bnuc.i + B nuc ,i)/2. Several ratios of dot nuclear fields, r — 
B nuc ,i/ B nuCtT are considered: r = 1 (black), r = 1/3 (red) 
and r — (blue). Cusp behavior near zero field is found in 
the limit of highly inhomogeneous dot field strengths. 



The mean values of the coefficients |c ++ 1 2 , |c_| | 2 are in 

turn st raig htforward to calculate approximately (as done 
in Ref. l2fj. supplemental information 3 ), giving 



r + _ = r( e )(| c+ _| 2 ) I 
r ++ = r( e )(| c++ | 2 ) I 



4 
T(e) 



(1 + IJr) (28) 
(1 - I^r) (29) 



Contrary to expectation, the blockade is contributed 
to solely by the \m s \ — 1 triplet states. In particular, 
spin blockade is charge transport at finite bias through, 
for example, the charge states (0,1) — > (1,1) — > (0,2). 
For biases between left and right leads that are less than 
J only the four (1,1) spin states and the state 1(0,2)5*) 
are necessary for understanding the process. An elec- 
tron (of arbitrary spin) loads from the left lead, creating 
with equal probability any of the states \s, s'). This then 
tunnels with a rate T + _ or T ++ to the state 1(0,2)5*), 
after which the extra electron on the right tunnels into 
the leads, and the cycle repeats anew. The average cur- 
rent through the device is dominated by the slowest rate, 
which in the absence of cotunneling, is In other 

words, loading into a spin-aligned state \s, s) prevents 
further charge transport until it decays, with rate 
or is replaced from the leads by a cotunneling process. 

The measurements by Johnson et al. demonstrate that 
the transition probability from (1,1) to | (0,2)5) depends 
strongly on both magnetic field and detuning 2 ^. Our 
theoretical model, which accounts for hyperfine mixing 
coupled with inelastic decay agrees well with experimen- 
tal results for timescales less than 1 ms. The discrepancy 
between experiment and theory for longer times suggests 
that other spin relaxation processes may become impor- 
tant above 1 ms (spin-orbit). 



with 



(30) 



and similarly for I r . We plot the product Iil r , found in 
both Eqns. l2"8ll29"l as a function of external magnetic field 
for increasing difference in dot sizes in Fig. [5] This indi- 
cates that the effective decay rates are largely indepen- 
dent of the ratio of dot sizes, relying only on the average 
effective nuclear field, B nuc = {B nnc ,l + B nuc ^ r )/2. 

We find that past the charge transition with e ^ \T C \ 
the states |5) and \Tq) decay to | (0,2)5) with a lifetime 
r^ 1 :, while the two \m s \ = 1 triplet states have a lifetime 

At finite magnetic field, T_| r ++ , and we can 

call the states of the |m s | = 1 subspace "metastable" . 
The metastability allows for charge-based measurement 
to distinguish between {|5),|T )} and \T±) subspaces: 
by using a nearby charge sensor, the decay of \S) , |T ) 
may be detected long before \T±) has finite probability 
of decay in the weak tunneling limit. This indicates that, 
while at zero field decay of the (1,1) states to the state 
|(0,2)5) is governed by a single exponential, a double- 
exponential behavior appears as B > B nuc is satisfied, in 
direct confirmation of the results of Ref. [20. 



3 In contrast to Ref. [201 . where B n 

Bnuc = \J (|5 nuc |2)/3. 



|Bnuc| 2 }, we define 



III. QUANTUM CONTROL OF TWO 
ELECTRON SPIN STATES 



We now analyze how time-dependent control of gate 
parameters (e.g., e) may be used to control electron spin 
in double quantum dots. Of particular interest are meth- 
ods for probing the hyperfine interaction more directly 
than in the previous section. The new techniques we 
use are primarily rapid adiabatic passage and slow adia- 
batic passage. Rapid adiabatic passage (RAP) can pre- 
pare a separated, two-spin entangled state (|5) in the far- 
detuned regime), and when reversed, allows a projective 
measurement that distinguishes the state |5) from the 
triplet states, \T m ). A similar technique used at large 
external magnetic field, slow adiabatic passage (SAP), 
instead prepares and measures eigenstates of the nuclear 
field, |s, — s). We connect these techniques with the ex- 
periments in Ref. [2l| and estimate their performance. 



A. Spin-to-charge conversion for preparation and 
measurement 

Adiabatic passage from e <C — T c to e ^> T c maps the 
far-detuned regime states |5) , \ T m ) to the states past the 
charge transition |(0, 2)5) , \T m ), allowing for a charge 
measurement to distinguish between these result o 21 ' 26 . 
In the quantum optics literature, when adiabatic trans- 
fer of states is fast with respect to the relevant dephasing 
(nuclear spin- induced mixing, in our case), it is called 
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FIG. 7: (a) Rapid-adiabatic passage: starting in the state 
|(0, 2)S), the detuning is changed from e 3> T c to e <C —2c, 
fast with respect to the nuclear energy scale, 7 e B n uc. (b) Slow 
adiabatic passage: as above, but once the system is past the 
S — T+ degeneracy point, the change of e is made slow with 
respect to the nuclear energy scale. The zoomed-in section 
shows the current nuclear energy splitting (ill) and the nuclear 
field eigenstates, \s, — s) and |— s, s). Both procedures may 
be reversed to transfer either |5) (RAP) or \s, — s) (SAP) to 
|(0, 2)S) while keeping the other states within the (1,1) charge 
configuration, allowing for charge-based measurement of the 
system. 



"rapid adiabatic passage" (RAP) and we adopt that ter- 
minology here. 

When the change of detuning, e, is adiabatic with re- 
spect to tunnel coupling, T c , but much faster than 7 e -B nU c 
(the hyperfine coupling), the adiabatic passage is inde- 
pendent of the nuclear dynamics. For example, starting 



past the charge transition with the state |(0, 2)5) ( 

and t » r c and using RAP to the far-detuned regime 
causes adiabatic following to the state IS*). This pre- 
pares a separated, entangled spin state. The procedure 
is shown in Fig. [7^,. 

The reverse procedure may be used to convert the 
singlet state to the charge state (0,2) while the triplet 
remains in (1,1). Then, charge measurement distin- 
guishes singlet versus triplet. Specifically, if we start 
with some superposition in the far-detuned regime, = 
°s \S) + J2m c m \T m ) where eg, c m are quantum ampli- 



tudes, after RAP past the charge transition, with e ^ T c , 
the state is \ f) = c s e l * | (0,2)5) + £ m c m \T m ), where 
4> is the adiabatic phase accumulated. Recalling that 
|(0, 2)5') is in the (0,2) charge subspace, while \T m ) states 
are in the (1,1) subspace, a nearby electrometer may dis- 
tinguish between these two results, performing a pro- 
jective measurement that leaves the \T m ) subspace un- 
touched. Furthermore this measurement is independent 
of the adiabatic phase. 

In contrast, if the change of e is slow with respect to 
nuclei, adiabatic passage follows to eigenstates of the hy- 
perfine interaction. For simplicity, we assume that RAP 
is used between the charge transition to just past the 
S — T + resonance of Eqn. [251 such that we may neglect 
transfer between the 5 and T + states (see also Fig. (SJd. 
This requires an external magnetic field |-B C xt| 3> -Bnuc- 
Continuing from this (e < — T c ) point to the far-detuned 
regime, e is changed slowly with respect to fi7 e -B n uc- Ac- 
cordingly, adiabatic passage proceeds into eigenstates of 
the nuclear field, |s, — s), as shown in Fig.[7)3. These are 
the product spin states, with one spin up and the other 
down with respect to magnetic field. This technique may 
be called slow adiabatic passage (SAP). 

Rapid adiabatic passage maps |5) <-> 



e** 1(0,2)5), \T„ 



\T m ), leaving the triplet states 



unperturbed. For SAP, the mapping is 



j<t> 



|(0,2)5) <- , 
\To) 

\T±i) <-> \T ±1 



<-> e**' \-s,s) 



(31) 



It is always the current, lowest energy eigenstate of the 
nuclear field that | (0,2)5) maps to. That is, we choose 
s such that E s .- S < -E-s.s, with E s .- S = sj e (B l I] 



B^ uc z ), (see Section |nl Eqn. [2H evaluated at large ex- 
We remark that SAP allows for 



ternal magnetic field). 

deterministic preparation and measurement of states ro- 
tated 7r/2 with respect to |5) , \Tq) without direct knowl- 
edge of which states they correspond to for each realiza- 
tion. 

We now examine the adiabaticity condition for slow 
adiabatic passage. In SAP, e is changed at a constant 
rate to approach the far-detuned regime (point S). Using 
the approximate relation J(e) = 4T c 2 /e, the adiabaticity 
condition is fi J <C (hw) 2 . Neglecting factors of order 

unity, this can be rewritten He^§- (hcu) 2 . 

As a specific case, we consider T c ~ 5 /ieV, and e G 
5)) [50,550] /^eV. The required time to make the 500 /ieV 



change, r, gives hi = ^500 ^eV, and roughly, he-^- = 
f x 50 neV. In units of time, H 2 (t x 13 ns). For r = 1000 
ns, the adiabaticity requirement is that the current value 
of 1/M < 3 x V1000 x 13 ns = 300 ns. For the nuclear 
fields in lateral quantum dots such as those of Ref. [20l 
(each dot with a T 2 * = 10 ns), the probability of l/|a->| > 
300 ns is 



P(\oj\ < 3 ^ 



3 /AS 



dv 



(32) 
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This gives an error probability of 3% for 300 ns rise time, 
that is, every 1 in 30 experimental runs, the nuclear gra- 
dient will be too small for the adiabatic filling of |±) to 
occur. 

We can ask the effect of finite, residual exchange en- 
ergy J rcs at the S point. Finite J leads to filling of a 
superposition of \s, — s) and |— s,s): 

cos(0) \s, — s) — sin(0) |— s, s) , (33) 

where the value of s is, as above, determined by the cur- 
rent value of lu and 4> — arctan[ J res / (-\/ J? es + u> 2 + \u>\)] is 
the adiabatic angle. The resulting loss of contrast will be 
sin 2 (</>) ~ (J/2uj) 2 . For residual J rcs ~ 0.1 fi, the error 
is less than 1%. For J rcs ~ 0.25 £1 the error is order 2%. 
The role of residual nuclear fields during the exchange 
gate is evaluated elsewhere^. 

B. Probing the nuclear field and exchange 
interactions with adiabatic passage 

We now consider how adiabatic passage can be used 
to probe the dephasing and exchange energy of a spin- 
singlet state. This relates directly to a critical question in 
quantum information science: how long can two electron 
spins remain entangled when the electrons are spatially 
separated on a GaAs chip. In our model, variations in the 
local nuclear environment cause the spatially separated 
electrons to experience distinct local magnetic fields, and 
hence precess at different rates, mixing the singlet and 
triplet states. If many measurements are taken to de- 
termine the probability of remaining a singlet, the time- 
ensemble-averaging leads to an observable dephasing of 
the singlet state (T 2 *)2i. 

To evaluate the effects of nuclei on this process, we will 
calculate the singlet autocorrelation function As{t) = 
| (S(t) \S(0)) | 2 for the far-detuned regime. This autocor- 
relation function has been evaluated for quantum chem- 
istry^, but not for this specific scenario. We remark that 
our approach is similar to the single dot case considered 
in Refs. [lOGJ. 

We start by evaluating the evolution operation U(t), 
where the Schrodinger picture \S(t)) = U(t)\S). Tak- 
ing J — > for the far-detuned regime, we solve ana- 
lytically the equation of motion any spin state of the 
(1,1) subspace. In particular, we take the Hamilto- 
nian of Eqn. [18] and write it in the form two effective 
fields, each acting separately on one spin. The evolu- 
tion operator, U(t) = exp(—iHt/K) can be factorized as 
U(t) = Ui{t) <g> U r (t), where 

Ui(t) = cxp(-i le t[B + S n uc,i] • Si) (34) 

is a rotation of spin i about an axis tii = (xi,yi,Zi) cx 
B + B nuCi i of angle tuji where 

UJi = J e \B + B nuC;i |/2. 

If the system is prepared by RAP in the state 
\S{t = 0)} = (HI) - ||T))/\/2, and subsequently mea- 
sured using RAP to distinguish the singlet and triplet 
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FIG. 8: Pulse sequences: detuning parameter e versus time, 
for (a) RAP for measurement of the singlet autocorrelation 
function, (b) SAP for the exchange-gate sequence, and (c) the 
singlet-triplet echo sequence. Blue is the charge transition 
region, while yellow is the far-detuned regime. The charge 
degeneracy point (the crossing from (1,1) to (0,2)) and the 

degeneracy between js 1 ^) and |T+) (when j(e) = 7 e |S ex t|) are 

shown for reference. 



subspaces, the measurement probes the state's autocor- 
relation function. Starting in a singlet at t = 0, the 
probability of remaining a singlet after a time t is given 
by the autocorrelation function 

A s (t,B) = \{S(t) \S(t = 0)) | 2 (35) 
= | cos(w/t) cos(u! r t) + fii ■ ri r sin(a;;£) sm(u> r t)\ 2 

To obtain the signal in the quasi-static approximation, 
Eqn. [35] must be averaged over the different possible nu- 
clear field values. We examine the zero-field and finite- 
field cases. 

When B = 0, the properties of tii within the QSA 
are described by (ni^rij^) — 5^5^/3. Averaging over 
nuclei, the signal is 

As(t,0) = (cos 2 (cj^)) nuc (cos 2 (w r i)) nuc + 

-(sm 2 (w/i)) (sin 2 (u) r t)) nuc (36) 
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where 

(cos 2 M)) nuc = l[l + e-^/ T l^(l-(t/T*.) 2 )], 
(sin 2 M)) nuc = l[l- e -^/ T ^ 2 (l-{t/T*.) 2 )]. 

We recall that = (7e-Bnuc,i) _1 - A distinct difference 
of this model from other dephasing mechanisms is the 
order 10% overshoot of the decay at short times, and 
the asymptoptic approach of As(t ^> T% , 0) to 1/3. A 
classical master equation would exhibit neither of these 
features — they are unique identifiers of the quasi-static 
regime, in which different coherent dynamics are aver- 
aged over many realizations. Numerically we find that 
these qualitative features do not depend on the relative 
size of the two quantum dots (B nuc j/ B nucr ) for varia- 
tions of up to 50%. 

Another regime of interest is when the external field 
is much larger than the effective nuclear fields (\B\ ^> 
B nuc .i)- Spin-flip terms are highly suppressed and the 
system is restricted to two levels, \S) and the m s = 
triplet, \T ) = (|U) + |4T))/V2. This is described by the 
Hamiltonian of Eqn.[2U This effective two-level system's 
evolution operator is straightforward to evaluate: 

A s (t,B > B nuc ) = (cos 2 (^)) = ~[1 + e ~(*/T 2 % ff ) 2 ] . 

(37) 

with T 2 * eff = l/^[(T 2 * i )- 2 + (T 2 * r )- 2 ]. Qualitatively, 

the decay of the autocorrelation function As due to 
the nuclear field is described by Gaussian decay with a 
timescale eff . Similar to the case of zero magnetic field, 
the behavior of this autocorrelation function is indepen- 
dent of variations in dot size up to ~ 50%. 

We now indicate how slow adiabatic passage at large 
external magnetic field allows measurement of the results 
of an exchange gate. In particular, SAP allows for prepa- 
ration and measurement of the individual spin eigen- 
states, |l/2,— 1/2) and | — 1/2,1/2). An exchange gate 
leads to partial rotation between these states, where the 
rotation angle is given by the product of the exchange 
energy during the gate, J(e), and the time the exchange 
energy is non-zero, Finally, reversing SAP takes the 
lower energy eigenstate (|s,— s)) back to 1(0,2)5*) while 
the higher energy eigenstate (|— s,s)) is mapped to \Tq), 
a (1,1) charge state. The final measurement determines 
whether the final state is the same as the initial state 
(the (0,2) result) or has changed to the state with the two 
spins exchanged (the (1,1) result). Thus, preparing the 
state 1 1/2, —1/2) and measuring in the same basis distin- 
guishes the results of the exchange- based rotation of the 
two-spin state. For example, when the probability is 50% 
for either measurement result, a \/SWAP gate has been 
performed. When the probability goes to 100% of recov- 
ering the higher energy eigenstate (measuring (1,1)), a 
complete swap of the two spins has occurred (SWAP). 

As before, we consider the probability of returning to 
the lower energy eigenstate. Now, however, this state is 



the |s, — s) state as described in the previous subsection. 
After preparing this state, we perform the resonant ex- 
change gate of angle Qe = Jt/H where t is the time spent 
waiting with exchange energy J. This leads to a rotation 
of the state |s, — s). Its autocorrelation function is given 
by 

A s ,_ s (t) = || s ,- s (t))( s ,- s || 2 (38) 
= cos 2 (9 E /2) . (39) 

If the exchange term J < fi7 e -B nuc , then additional effects 
due to nuclei would be observed; we evaluate these below. 

We emphasize that the combination of RAP for prepa- 
ration (prepares \S)), SAP for preparation (prepares 
lower energy eigenstate of the nuclear field, |±)), RAP 
for measurement (spin-to-charge in \S) , \Tq) basis) and 
SAP for measurement (spin-to-charge in current eigen- 
states of the nuclear field, |±)), when combined with the 
exchange gate (rotations of |±) to |±) + i |=p)) allows for 
full state tomography in the IS*) , \Tq) subspace. 

C. Errors in exchange gates 

The primary error in exchange gates is likely due to 
charge-based dephasing and is directly related to the 
parametric dependence of the exchange energy, J, on gate 
voltages near the charge transitio n 27 i 30 i 55 . In addition, 
other errors are possible due to the stochastic nature of 
the nuclear field. For example, there is the possibility 
that the current value of the field gradient frequency, u>, 
is sufficiently small to make the initial and final transfer 
stages non-adiabatic. Also, the gradient can flip sign in 
the course of the experiment. Finally, finite residual ex- 
change interaction during SAP reduces effectiveness. We 
consider each of these in turn below. 

In the far-detuned regime, the energy gap between the 
(1,1) singlet/triplet space and higher orbital states, as 
well as (2,0) and (0,2) charge states, is large (of order 
hujo, the orbital level spacing of a single dot, and E c /2, 
the single dot charging energy, respectively). At dilution 
refrigerator temperatures, this energy gap is many times 
greater than k^T, and absorption of a quanta of energy 
from the environment leading to incoherent excitation 
may be neglected. Also in this regime, the residual ex- 
change splitting is both numerically small and insensitive 
to first order fluctuations in detuning, e, leading to little 
charge-based dephasing due to differential coupling of the 
\S) state to the doubly-occupied states when compared 
to the \T m ) states' couplings to doubly-occupied states. 
The system remains sensitive to charge-based dephasing 
up to second order due variations in the tunnel coupling, 
T c . If we can write variations of T c from the mean as 8T Cl 
its correlation function is given generally by 

(5T c (u + t)5T c {u)) =h 2 J dw S Tc (u)e iuT • (40) 

The corresponding phase noise term in the Hamilto- 
nian is V n = vST c (u)(J2 rn \T m )(T m \ - |5)(5|)/2 where 
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(a) 



n8/ 



- |T-> (b) |T-> ■ 



-|G> 



-|T > 
IS) 



|T+) 



|To> 

r)'8f+ri"6e- 
|T+> 



r(i) 



is) 



Far-detuned 



Biased 



FIG. 9: Level structures of the double dot system in the 

(a) far-detuned regime, with charge dephasing shown, and 

(b) near the charge transition with charge-based decay and 
dephasing. 



5 meV, J = 0.3 



again using E c 

meV, and (T c ) = 0.01 meV. 

As an example, a coherence between the \S) and \T m ) 
subspace in the far-detuned regime could be expected to 
decay due to the noise on T c according to 

( 41 ) 

While we consider a variety of noise sources in Ap- 
pendix [B} it is instructive here to take the case of white 
noise spectrum with St c — 7o/27r. It leads to exponen- 
tial decay of coherences between IS*) and states of the 
\T m ) subspace with a constant 7t c = rj 2 ^. In general, 
as i] <§; 1, this decay will be negligible. 

The charge transition will have stronger dephasing 
when compared with the far-detuned regime. In addi- 
tion to inelastic decay of the excited adiabatic state to 
the ground adiabatic state, the system is susceptible to 
fluctuations in both e and T c , as J is potentially large 
and strongly dependent on gate parameters. In so far 
as the power spectra associated with T c and e have no 
appreciable spectral components at frequencies of order 

G) and 



e, the excitation/relaxation terms between Gj and S 

may be neglected. Then we may restrict considerations 
to dephasing of coherence between the subspace of \T m ) 

and the state S ) , and write the effective Hamiltonian as 



Vu-02 = [v'ST c (T)+ V "Se(T)][y2 \T m )(T m \ - \S)(S\}/2 



with rf = 



2T r 



(42) 

and rj" = |(1 + 7J ^ l ). The ac- 
companying noise is given again by Eqn. 1411 but with 
f] = \J (r/) 2 + (j]") 2 and the power spectrum replaced 
with 

Sii_o2M=(~) StM+^j) (43) 

While rf may be small (rj <C 1), near the charge tran- 
sition rf' is of order unity. This indicates that noise in 
the detuning parameter e has significant repercussions 
for coherences between singlet and triplet states during 
exchange gate operation. 




10 

t[ns] 



15 



20 



FIG. 10: Decay of exchange oscillations for the four scenar- 
ios (clockwise from lower left): white noise, 1//, ohmic, and 
super-ohmic (co 2 ). White corresponds to probability 1 of end- 
ing in the initial state \s, —s) after exchange interaction is on 
for a time te (bottom axis), while black is probability of 
ending in the initial state. Tunnel coupling is taken to be 10 
/xeV, and the spectral density (y) is chosen for similar behav- 
ior near e = 50 /ieV. Note that 1/f terms increase decay in 
the slow oscillation limit, while increasing powers of u) (white 
noise=0, ohmic = 1, super-ohmic = 2) lead to more oscilla- 
tions for smaller exchange energies. 



We consider charge-based dephasing for the exchange 
gate according to Eqn. SU In all cases we assume the 
spectral function has a high-frequency cutoff 7 such that 
7 < J. This assumption prevents dephasing noise from 
producing population changes (relaxation) due to en- 
ergy conservation. Additionally, we rewrite the expected 
probability of maintaining phase coherence (Eqn. I4ip as 



exp(— T] 2 P) 



(44) 



where P = J du; S{u) a ™\^ 2 "J 2) is set by the spectral 
function and the time of the exchange gate, te- This 
allows separation of the interaction strength (rf) and the 
noise spectrum. Each spectrum considered has a normal- 
ization parameter v such that S(u>) has units of inverse 
time. 

We note that in general, the number of observable ex- 
change oscillations will be limited by these dephasing 
processes. By finding T2 = P _1 (t] _2 ) the observable 
number of oscillations goes as TzJ/h. When, for exam- 
ple, S(u>) — r/27r, we may easily invert P(t) and find 

t 2 = v - 2 /t. 

Comparing the behavior of the ohmic and super-ohmic 
cases to the 1/f case (see Appendix [Bjl . the limiting value 
of P for the super-ohmic case and the power law tail of 
the ohmic decay indicates that for small coupling pa- 
rameter 77, the super-exponential 1/f terms (going as a 
gaussian) will dominate at long times. For very short in- 
teraction times, all three will be less than the white noise 
contribution, which goes linearly in time. The different 
behaviors are shown in Fig. [TUJ 

This indicates that electrical control of exchange inter- 
actions in double dot systems may be relatively robust 
with respect to nuclear spin degrees of freedom. How- 
ever, during the exchange gate, the system is susceptible 
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to charge-related dephasing. The observed decay of os- 
cillations of Ref. Hi], in which the decay rate appears to 
match the exchange energy such that the number of ex- 
change oscillations observed is independent of the detun- 
ing, is qualitatively similar to the behavior of oscillations 
in the presence of sub-ohmic noise. A more detailed ex- 
perimental analysis of the noise will be required before 
a direct comparison between experiment and theory will 
be possible. 



IV. EXCHANGE GATES AND ECHO 
TECHNIQUES 

The techniques of rapid adiabatic passage, slow adia- 
batic passage, and time-resolved control of the exchange 
interaction in the previous section have far reaching ap- 
plications for quantum control of spins^S. In this section, 
we consider how exchange gates can undo the effect of the 
quasi-static nuclear field, greatly reducing the deleterious 
effects of nuclear spins on electron spin coherences and 
allowing for electron spin measurements to determine nu- 
clear spin correlation functions. 



A. Spin-echo in the singlet-triplet basis 

Since the nuclear fields vary slowly on timescales com- 
pared to a typical pulse cycle time, a spin-echo pulse 
sequence can be used to refocus the spin singlet state. A 
spin-echo pulse sequence based on fast electrical control 
of the exchange interaction was demonstrated in Ref. [2l]. 
The experiment starts by using RAP to transfer |(0, 2)5') 
to 1 5), preparing a separated singlet pair of electron 
spins. As demonstrated in the T£ experiment, the hy- 
pcrfine interaction mixes the singlet and triplet states 
on a 10 ns time scale. After a separation time ts, a it 
exchange pulse (SWAP) is performed by adjusting the 
detuning to a region with a finite exchange energy. The 
exchange energy is then set to zero by moving to the far- 
detuned regime for a time ts>, during which the singlet 
state refocuses. 

The dephasing due to hyperfine interactions occurs by 
producing a relative, unknown phase between |||) and 
lit)- Switching between these two states via exchange 
gates will produce an echo (recovery of the original state, 
\S)) if the waiting time before and after the SWAP op- 
eration is the same. In other words, the pseudo-spin of 
IS*) and |T ) (the m s — subspace) is amenable to echo 
techniques using exchange interactions. 

We will use the pulse sequence of Fig. [3fc. In the far- 
detuned regime at finite field, 15*) and \Tq) mix due to 
nuclei, as per the Hamiltonian of Eqn. [5T] This mixing 
is driven by a relative constant, but unknown, energy, 
corresponding to the current value of 



Within the corrections to the quasi-static approximation, 
lu = jeiSBi.zit) — 5B rtZ (t), where 5B^ r ^ z {t) are random, 
Gaussian variables described by Eqn. [T^l For clarity, we 
rewrite it here: 



(v(t + r)w(t)) = duj S(uj) 



(46) 



For example, the singlet-singlet correlation function at 
finite field (Eqn. |3~7| is modified, noting that 



(cos 2 (/ uj{t')dt')) 
Jo 



1 + exp(— J duj S(oj) 



sin 2 (™/2) 



(47) 

For S{lo) with a high frequency cutoff below 1/r, we 
can Taylor expand the sin term. Then, comparison with 
Eqn. [37] indicates: 



(48) 




Now, we consider how this result changes with the 
more complex sequence given by Fig. [5J:, in particular 
its dependence on te and Tgi . If te = 0, then nothing 
has changed from before, except now the system evolves 
according to U nuc (T S >)U nuc (T S ) = U quc (t s +t s >). How- 
ever, what happens with finite r^;? 

For our reduced two-level system, when te 7^ the 
effective Hamiltonian during this stage is given by 



H e ff,E = hwa x /2 + J{tE)<J z /2 



(49) 



hw = 7 e (B, 



IlUC,/,2 



B n 







(45) 



where Pauli matrices with states \S) , \Tq) as a psuedo- 
spin, in our logical basis defined above (i.e., a z = 
|To)(T | - \S)(S\). Taking the turn-on of finite J to be 
instantaneous, the total evolution operator is 

U(t S ',te,ts) = U auc (T S ')exp(-iH e ff iE T E /h)U BUC (T S ) . 

(50) 

When the exchange energy for the middle point satisfies 
J{ce) 7eB nuc , we approximate the middle term of the 
evolution operator by Ue{te) ~ exp(— ite J(eE)cr z /2h). 

For now, taking Co to be constant, at the end of the 
pulse sequence, the probability of returning to the state 
|(0,2)5) (C) is given by 



C = Tr nuc (0| e -^s>zaj2 e -iT E j(e E )az/2h e -iT S zaj2 ^ 

(51) 

To see the analogy between this evolution and a spin echo 
experiment, we insert unity after the initial ket, i.e., |0) = 

e iJ(e E )a z /2h e -tJ(e E )(-l)/2h |q^ &g q ^ ^ _j eigenstate f 

o z . The overall phase is irrelevant due to the absolute 
value terms, and so we rewrite the above as 
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C = Tr„ 



(52) 



The term in parenthesis, an exchange gate, is a rotation 
of angle 9e — teJ(ze)/Ti about the z axis of the psuedo- 



spin. It acts on the operator U nuc (Ts) = e 



.flip- 



ping the sign of the a x operator when 9 E = n (SWAP). 

We can probe its effects by analogy with standard spin 
echo. For example, when 9e = 7T, we get C = ^[1 + 
cxp(— (u) 2 )(ts — Tg') 2 /2)]. This means the dephasing due 
to nuclei in the first half of the sequence is exactly undone 
when the rephasing time is equal (rg = rg/). For fixed 
T s + T S' j the probability of returning to the singlet state 
as a function of t$— ts> and te should exhibit 50% mixing 
expected for dephasing when \rs — ts>\ > T 2 * , and when 
teJ '/h is not an odd integer multiple of n. For example, 
setting ts + ts' — 100 ns, and using ce such that J{€e) = 
/ry e (20 mT), the probability exhibits this behavior. 

In practice, the instantaneous approximation breaks 
down in realistic situations, as does the J 3> 7 e -B n uc as- 
sumption. The former is easily fixed by noting that such 
a "Rabi" type pulse only has sensitivity to the integrated 
area, i.e., 9 = J J(e(t))dt. The latter requires working 
with finite r^fL This effect has been considered in detail 
elsewhere^. 



B. Probing nuclear spin dynamics through echo 
techniques 

So far the analysis has worked entirely within the QSA. 
However, the echo sequence in principle reveals higher 
order dynamical information about the nuclear field and 
other noise sources. Effects of nuclear spin dynamics on 
electron spin decoherence has been considered by sev- 
eral author o 14 i 15 ' 17 i 28 ' 29 . We now consider the echo se- 
quence with to(t) a Gaussian variable slowly varying in 
time. This allows a large range of possible noise sources 
and correlation functions to be considered. 

Assuming (for the moment) that the exchange gate is 
of high fidelity and insensitive to the current value of A z , 
the measured result of the echo sequence is given by 



rt s +t s , 

C = |(0 L |exp(z / u)(t)dta x /2 
Jt s 

[•ts 

-i / to{t)dta x /2) |0 L ) I 2 



(53) 



ts+t s , 



cos 2 [( uj(t)dt+ Cj{t)dt)/2] (54) 

Jt s JO 

1, 



|e - a + e l = + 2] 



(55) 
(56) 



where S = J. 
mcnt of S is 



t s +t t 



' Co{t)dt + J * s Co{t)dt. The second 



mo- 



/~2\ 



(si, 



(3D 

4 „• 2/ 



2(H<j,H 



(57) 



with (S 2 ,) = Jdu S(v)J I sm 2 (t s ,v/2) ~ (t s < /T 2 *) 2 and 
similarly for S5. The cross term, corresponding to the 
correlations between the two frequencies, is 

2(Z S ,Z S )= I dv g(,) Sinfeiy ( 2) s \f W2) e » (ts+ts , W2 

J W 2 ) 

Finally, when t$ = tgr, the second moment is 
sin 2 (t s ^/2) 



= 2 dv S{v) 



1 



Jt s v\ 



(59) 



W2) 2 

For low frequency noise, with cutoff 7 ~ 
[jdu ^Siv)] 1 ' 2 < 1/ts, we obtain (S 2 ) ~ 
2(ts) 2 x (ts'y) 2 , or decay in the total wait time 
2ts with an effective time constant 

3 V4 (T * /7) V2 



T, 



2.SE 



(60) 



Our evaluation has implicitly assumed that above the 
cutoff, 7, S(uj) dies off at least as 1/w 3 . We suggest that 
this is appropriate for interaction times on the order of 
microseconds — initial decay of the nuclear spin correla- 
tion function is quadratic. For longer interaction times, 
a different decay morphology (going as exp(— t 3 )) could 
be observed. 

Our predicted deca y is consistent with the experiments 

of Ref.dH. In Refs.^jlMHIil 7 is of order 10 1118-1 ' S iv " 
ing T2,se = 2 [is. However, addition exponential decay 
could be observed if the high frequency cut-off assumed 
above has a too-slow decay, going as cj -2 ~ £ with e < 1. 
Furthermore, higher order pulse sequences, such as Carr- 
Purcell, will likely allow for extensions of the echo signal 
to substantially longer time a 57 i 58 ' 59 . 

The experiments discussed in the previous sections 
demonstrate that the hyperfine interaction is efficient at 
dephasing an initially prepared spin singlet state on a 10 
ns timescale. By using a simple spin-echo pulse sequence, 
this bare dephasing time was extended by over a factor 
of 100, to times T2>1.2/is. Further experimental effort 
will be required to fully map out the nuclear correlation 
function and extend the lower bound on electron spin 
coherence times. 



V. CONCLUSIONS 

We have shown how a model combining charge and 
spin interactions for two-electrons in a double quan- 
tum dot effectively describes the experimental results 



15 



of Refs. I20II21I . By starting with the case of a single 
electron in a single dot, we employed the quasi-static 
approximatio n 10 ' 12 ' 15 ' 17 ' 18 ' 27 ' 28 ' 29 , and considered its ap- 
plicability to describing current experimental results. 

The spin interactions with nuclear spins are extended 
to the double dot case, and two regimes emerge: the far- 
detuned regime, in which two electrons are in separate 
dots and interact with independent nuclear fields, and 
the charge transition, in which the two electrons may 
transition from a separated orbital state to a doubly- 
occupied, single-dot state. 

This model was used to describe spin blockade, and we 
found spin blockade is broken by interactions with nuclei 
near zero magnetic field, explaining the experimental re- 
sults of Johnson et al. (Ref. I20I ). A striking magnetic 
field dependance is derived, consistent with observed ex- 
perimental behavior. This indicates that the dominant 
mechanism for spin blockade at finite magnetic field is 
trapping of the m s = ±1 separated spin states, as their 
mixing with the charge-transition-allowed singlet \G) is 
suppressed by finite Zeeman splitting. Observation of the 
breaking of spin blockade near zero field provides a sen- 
sitive measure of the magnitude of the random nuclear- 
spin-induced magnetic field. 

Time-domain control of local potentials, achieved 
through manipulation of electrostatic depletion gates, 
provides powerful mechanisms for preparing and mea- 
suring spin singlets, as well as eigenstates of the nuclear 
field. These techniques have been exploited by Petta et 
al. to measure the effective dephasing of a two-spin en- 
tangled state, and to probe via coherent oscillations in 
the m s = two-spin subspace the exchange interaction 
as controlled by gate voltages^ 1 . Limiting mechanisms 
for such oscillations, due to charge fluctuations of inde- 



terminate nature, are considered for a variety of environ- 
mental noise spectra. 

Finally, we analyzed how controlled exchange interac- 
tions can protect the electron spin from the deleterious 
effect of nuclear spins by working within a two-spin sub- 
space, putting in specific terms protocols previously con- 
ceived more generally. We expect the limiting mechanism 
for the rephasing of the two-spin states comes from cor- 
rections to the quasi-static approximation, and as such, 
spin-echo experiments provide a useful measure of the 
validity of this approximation. 
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APPENDIX A: ADIABATIC ELIMINATION FOR 
NUCLEAR-SPIN-MEDIATED INELASTIC 
DECAY 

We will transform the superoperator (Eqn. |23|) into the 
interaction picture, but first introduce matrix elements 
between the eigenstates of the quasi-static fields and the 
state 15) occurring in the superoperator. For a single 
spin in a magnetic field B — B nuc (x,y, z) and \B\ = 
B nuc n (the roman variables x, y, z, n are chosen such that 
the nuclear field contribution will be of order unity), the 
eigenstates may be written by rotation from spin states 
aligned with the z-axis ({|t) , |l)}) : 



|l/2) 
-1/2) 



lim 



1 ( n + z z'-M/W|i) 

yj2n{n + z) \x' + iy -n - z I \ M ) 



(Al) 



The limit is taken only to remove the degenerate case of field anti-aligned with the z-axis, i.e., B = (0, 0, —B), which 
would be degenerate for this matrix, and is implicit in what follows. The corresponding eigenvalues of the Hamiltonian 
are ±Mln/2 in this notation. 

Setting \S) = (|Ti) — |4t))/"v2 and using the single spin transformations for I and r separately (with B l = 
B l nuc (xi,yi, zi),\B \ = m and similarly for r), we write 



c l/2,l/2 
c l/2,-l/2 = 
C-l/2,1/2 = 
c -l/2,-l/2 = 



1 1 

2' 2 
I 1 

2'~2 
1 1 
~2' 2 
1 1 

2'~2 

and Af 



1 

N 
1 



[{xi + iyi)(n r + z r ) - [m + z t )(xi + iy{)) 
[-(ni + zi)(n r + Zr) - (xi + iyi)(x r - iy r )] 

'1/2,-1/2 



c l/2,l/2 

y/8nin r (ni + zi)(n r + z r ) 



(A2) 

(A3) 

(A4) 

(A5) 
(A6) 
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It is convenient to define c ++ = 01/2,1/2 an d = C1/2,— 1/2 as the spin-aligned and spin-anti-aligned coefficients, 

respectively. This allows us to express \S) occurring in Eqn.[23]in terms of the eigenstates of the Hamiltonian as 



S,S = Ca 



1 1 

2' 2 



+ c+- 



1 1 

2' _ 2 



1 1 

2' 2 



1 1 

2'~2 



In the interaction picture, each eigenstate \s, s') evolves according to 



At)) 



— p -iE a .rt/h I J\ — „-'i(sn ! 7cS nuc , i +s'ri r 7 c B nuCir )t 



a nd\S(t))=J2 S!Sl c SlS ,\s )S >(t)) : _ 

With these results, the louivillian may be put into the interaction picture: 

P = T(e)/2 c.sZ,^-'*"' 3 ™^-^™' 1 ^'* ]t L s ,s',rAp\ 

s,s' ,r,r' 

where L.^ffl = [\s, s')(r, r'\p + p\s, s')(r, r'\ - 2\G)(r, r'\p\s, s')(G\] 



(A7) 

(A8) 

(A9) 
(A10) 



So far, this result is exact within the QSA. 

When r(e) <C "f e B nnc ,h JeB nuc ,r, the exponential phase terms of Eqn. lA9l oscillate substantially faster than p evolves. 
Adiabatic elimination becomes an appropriate approximation when we may neglect quickly rotating terms, i.e., if we 
may neglect certain degenerate cases, such as situations with \nij e B nuCt i — n r j e B nuc ^ r \ < T(e). In addition, we are 
implicitly assuming that r(e) ~ T(e + E SjS > — E r _ r i), which is appropriate for large e and smooth phonon density of 
states. 

More explicitly, we can average each term of over several spin rotations and make a Born approximation: 



,.[(r- S )»,7 s B m c ! ,+( r '- S ')» r 7 e B..c,r]t isjS , r/ ^ f )] 



2 i[(r-s)n i 7 e Bnuc,(+(»-'-s')™r7eB ml c, r .]t' c ^' 



The time averaging for a given s, s' , r, r' is straightforward so long as niB nuc i ^ n r B nuCt7 . 4 , giving 



lim — 

T — >oo T 



Jl{r-8)nn e B ano , l+ (r'-s')n r y e B Dnc , r ]t' dt > = ^ ^ ^ 



L s ,s' ,r,r'[p~{t)] 

(All) 
(A12) 



Thus terms with quickly varying phase go to zero. 



APPENDIX B: DEPHASING POWER SPECTRA 

We now evaluate dephasing in exchange gates due to 
charge fluctuations for a variety of spectral functions. 
The error should go as [1 — exp(— ^ 2 P)]/2, where the 
value 77 depends only on the detuning parameter. The 
impact of the particular spectral function is encompassed 

in P = J duj S(uj) S '"J/^ 2 ^ ■ We have assumed that S(lj) 
has a high frequency cutoff 7 <C \jt. 

1. White noise 

We set S(uj) — ^e^"/ 7 . Then we can evaluate 

P = f / du, *p*Me~h « vt . (Bl) 

This indicates the expected exponential decay of coher- 
ence due to white noise dephasing. 

4 When integrating over nuclear spin degrees of freedom, this 
corresponds to a surface of measure 0. 



2. 1// noise 

With S(u>) = v 2 /oj and frequency cutoffs B < uj < 
7 « 1/t, 

P = 21og[ 7 /P](i4) 2 . (B2) 

For a bath of 1// distributed fluctuators, the ini- 
tial dephasing is quadratic in the time of interaction, 
and increases as the measurement timescale (1/P) in- 
creases. At long times, the decay is gaussian, with super- 
exponential suppression of coherence. 



3. Ohmic noise 

Taking S(u>) = gtue^^^ 1 , evaluation of P is possible, 
giving 

P = 2 5 log[l + ( 7 t) 2 ] . (B3) 
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When considered in the decay function exp(— r/ 2 P), this 
gives a non-exponential decay law, exp(— r] 2 P) = [1 + 
(-yt) 2 ]~ 29n2 . In the short time limit, this is quadratic 
decay, going as 1 — 2(g , q 2 ){^t) 2 + O^-jt) 4 ), while the long 
time behavior is a power law with power —Agr/ 2 . 



4. Super-ohmic noise 



For the final spectral function considered here, we set 
S(u>) — v^^w^e - ^' 1 where £ > 1 indicates super-ohmic 
noise. Evaluation of P proceeds in a straightforward 
manner, giving 



P = 



i-C 



{' 



l + ( 7 t)3-§ cos[(C- l)tan" 1 ( 7 i)]| .(B4) 



where r(x) is the gamma function. This type of decay 
has a limiting value of 



lim P = 4T(C - 1) 

t^oo yy 

and short time behavior according to 



(B5) 



p(*)=2C(C-i)r(c-i) 



i-C 



( 7 f) 2 +0(( 7 t) 4 ) . (B6) 



For visual comparison, we calculate the expected, ob- 
servable Rabi oscillations using SAP as a function of time 
at finite exchange (tg) and at detuning e in Fig. [TU] 
In essence, increasing the exponent of the noise spectra 
(from 1/ui to constant to lu 1 ') leads to more oscillations 
as detuning is made more negative, i.e., as the admixture 
of charge decreases. 
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